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In this paper we study the jet response (particularly azimuthal anisotropy) as a hard probe of the 
harmonic fluctuations in the initial condition of central heavy ion collisions. By implementing the 
fluctuations via cumulant expansion for various harmonics quantifled by e„ and using the geometric 
model for jet energy loss, we compute the response Xn = Vn/en. Combining these results with the 
known hydrodynamic response of the bulk matter expansion in the literature, we show that the 
hard-soft azimuthal correlation arising from their respective responses to the common geometric 
fluctuations reveals a robust and narrow near-side peak which may be the explanation of the "hard- 
ridge" observed in experimental data. 
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I. INTRODUCTION 

The structure and properties of QCD matter under 
extremely hot and/or dense conditions are of fundamen- 
tal interest and provide unique environments for study- 
ing the strongest force of Nature. The hot deconfined 
QCD matter, the so-called quark-gluon plasma (QGP), 
was part of the history for cosmic evolution after the Big 
Bang and has now been created via relativistic heavy 
ion collisions (the "Little Bang" ) and explored in labora- 
tory experiments at the Relativistic Heavy Ion Collider 
(RHIC) and the Large Hadron Collider (LHC). In such 
collisions, highly energetic jets born from initial hard 
scattering provide natural "tomography" of the created 
hot QCD matter. Jet quenching due to energy loss along 
the jet path through the medium encodes essential in- 
formation about the dynamics of jet-medium interaction 
and the medium properties as well, which shall be in- 
ferrable from experimental observables such as the high- 
Pt hadron suppression and azimuthal anisotropy, as well 
as the di-hadron correlations (for reviews see e.g. [l|). 

While the jet quenching has been experimentally estab- 
lished as a very robust phenomenon at RHIC and now at 
LHC, the microscopic mechanism of jet energy loss is not 
yet fully understood. A conventional observable for quan- 
tifying jet quenching is the nuclear modification factor 
Raa which compares the particle production in the AA 
collision to the naive expectation from simply scaling up 
single NN cross section by the binary NN collision num- 
ber. A measured Raa significantly smaller than unity for 
the high-pt hadrons implies strong in-medium energy loss 
for the jets: indeed for central collisions at both RHIC 
and LHC we've seen Raa < 0.2. The Raa provides di- 
rect information on the average opaqueness of the created 
hot medium, and it is customary in various jet quench- 
ing models to use Raa in the most central collisions to 
normalize their respective parameters for the average jet- 
medium interaction strength. More sensitive are the geo- 
metric features of jet quenching observables that are par- 



ticularly useful in discriminating different models of en- 
ergy loss. These include the ^-dependence (when chang- 
ing colliding systems), the 5-dependence (when chang- 
ing the collision impact parameter or centrality class), 
and the 0-dependence (when changing the probe jet's az- 
imuthal orientation with respect to the reaction plane). 
For any given model with its parameters fixed in the most 
central collisions, the above geometric dependence and 
the correlations among different observables then provide 
crucial tests of the model P-[l0j. 

Let's elaborate a bit on the (/)-dependence of Raa- In 
non-central collisions, the medium "thickness" as seen 
by a penetrating jet depends on the azimuthal angle <f> of 
the jet with respect to the reaction plane, therefore lead- 
ing to the reaction-plane dependence of high-pt hadron 
suppression i.e. Raa{4>) 0- The dominant anisotropy 
in Raa{<P) (for non-central collisions) can be attributed 
to the second harmonic term cos(20 — 2"^ n p) with its 
coefficient being the elliptic "flow" parameter for high- 
Pt hadrons, V^"'^'^ which is a non-collective component 
though [11]. Despite the success of many models in de- 
scribing Raa and its centrality dependence, most mod- 
els significantly under-predicted the Vj''"'^'^ and failed the 
test by geometric data 3] 9]. The lack of a simultane- 
ous description for Raa and V^"'^'^ in a single model was 
not resolved till a new insight suggested in m. Moti- 
vated by the "magnetic scenario" for sQGP [ij], the au- 
thors of [3] pointed out that the energy loss of a jet may 
not simply scale with the local medium density as most 
models have assumed, but actually has nontrivial depen- 
dence on matter density (or temperature). It was par- 
ticularly shown that the geometric data Raa and ¥2°'^'^ 
versus centrality can be successfully described together 
by including a jet quenching component with strong en- 
hancement in the near-Tc matter by a factor 3^5 com- 
pared with higher-T QGP. Such an enhancement of jet- 
medium interaction may originate from non-perturbative 
structures created by the (color-) electric jet passing a 
plasma of (color-)magnetic monopoles that dominate the 
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near- Tc matter 12|[l^. A natural prediction of this sce- 
nario is that the effective jet-medium interaction would 
be rapidly reduced when going from RHIC to LHC en- 
ergies, or the hot matter created at LHC will be more 
"transparent" (apart from the trivial density factor) to a 
penetrating jet at LHC as compared with that at RHIC. 
Interestingly, recent quantitative computations and com- 
parison with LHC data at = 2.76TeV [14] 15] [8] in- 
deed suggest that 1) fixing the jet-medium interaction 
at RHIC and applying it directly to LHC would lead 
to "over- quenching" as compared with data; 2) reducing 
the jet-medium interaction by about a factor of 2 would 
allow a good description of the LHC data — such reduc- 
tion is remarkably rapid provided only a 30% increase in 
temperature from RHIC (^s = 0.2TeV) to current LHC 
{y/s = 2.76TeV)! The future LHC heavy ion data at 
y/s = 5.5TeV will be essential in establishing how much 
and how fast the jet-medium interaction decreases with 
temperature. 

In this paper we explore a new and interesting geo- 
metric aspect of jet quenching: the hard probe of the 
geometric fluctuations (in terms of various harmonics in 
aziniuthal angle) in the initial condition of heavy ion col- 
lisions. It has recently been shown in measurements [l6l - 
[20j and demonstrated in various modelings [2l| - l32| that 
there are very strong fluctuations in the initial matter 
profile from event to event. Such fluctuations contain var- 
ious higher order harmonics in azimuthal angle ~ cos(n(/)) 
rather than the naive expectation of 2n(i-harmonic dom- 
inance in averaged geometry. In the collective expansion 
of the bulk matter, these fluctuations lead to observed 
harmonic flows up to about n — 6 and also explain the 
soft di-hadron azimuthal correlations (the "soft- ridge" ) . 
A natural question to ask is, therefore, how a penetrating 
jet responds to the strongly fluctuating matter density 
from event to event. There have been some discussions in 
the literature [1] , but a systematic and clear picture 
has not been achieved. Intuitively one would expect the 
event-by-event azimuthal distribution of high-pt hadrons 
shall also reflect such harmonic fluctuations in the initial 
condition and it is our purpose to systematically quan- 
tify such response. In addition, both the quenching of 
hard jet and the expansion of soft matter are commonly 
correlated to the same underlying particular matter pro- 
file in a given event and such hard-soft correlation will 
survive the event average. We will examine how such 
correlation contributes to the di-hadron azimuthal corre- 
lation with a hard trigger and a softer associate hadron, 
in p articular a possible explanation of the "hard-ridge" 
[34| for which the origin has so far not been fully under- 
stood [11] [3^ [33. To highlight the role of pure fluctu- 
ations in leading to azimuthal anisotropy, we limit our 
study in the present paper only for the most central col- 
lisions, and to be specific we do calculations for RHIC 
i/s = 200GeV collisions. Unlike in non-central collisions 
with strong anisotropy already in the average geometry 
and dominated by the 2nd harmonics, in the perfectly 
central collisions (6 = 0) the average background geome- 



try is isotropic and the anisotropy from various harmonic 
fluctuations will be best manifested. 

The paper is organized as follows. In section II we 
will describe our model setup for the present study. The 
simulation results for jet response to the harmonic fluc- 
tuations in central collisions, characterized by Vn/f-n for 
n = 1,2,..., 6, will be presented in section HI. Such re- 
sults will be used in section IV to compute the hard- 
soft di-hadron correlations due to harmonic fluctuations 
where we will find a robust and narrow near-side peak in 
the azimuthal angle A(f> dependence as a possible expla- 
nation of the so-called "hard-ridge" . The summary and 
discussions will be given in section V. 



II. MODEL SETUP 

In this section we introduce our model setup. First we 
will show how we implement various harmonic fluctua- 
tions on top of the isotropic flreball in a perfectly central 
collision by using a modifled cumulant expansion based 
on that of 



25[. In particular we will discuss the fluc- 



tuations in both the participant density and the binary 
collision density profiles and their mutual relations. Sec- 
ond we will briefly discuss the geometric model we use 
for calculating the jet energy loss and its azimuthal angle 
dependence for a given matter distribution. 



A. Parametrization of Harmonic Fluctuations in 
Central Collisions 

In each central collision event, the participant density 
{pp) and the binary collision density {pc) are related to 
the thickness functions of the two colliding nuclei {Ta &): 
Pc(r) = aTair)Tbir), Pp{r) = Ta{r) [1 - Pb(r)] + (a 06). 
Here a{= 42 mb) is the N-N inelastic scattering cross 
section at y/s = 200GeV, A is the number of nucleons 
in each nucleus, and Pf, = (1 — aTb/A)^. It should be 
emphasized that in general one shall use Ta{f— b/2) and 
Ti,{r + b/2), and the above is only to be used for the 
perfectly central collisions with 6 = 0. Suppose event by 
event, Ta.t fluctuate around an isotropic background Tq, 
i.e. Ta,6 = Pq b + STa,b- As the result, by using the linear 
approximation, the fluctuations of pc and pp are directly 
related, as shown by 
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(1) 



(2) 



T° is obtained by the (optical) Glauber model. To pa- 
rameterize the fluctuations, we focus on Spp and later use 
the above equations to compute the corresponding Spc- 

We will use a modifled cumulant expansion based on 
the methods in [l^Hsj. Generally, we can make a Fourier 
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analysis for any profile 

Pp(r) = Pp,o{r) + 2pp,„(r) cos - (3) 



The difficulty is to parameterize the function form 
for the coefficients Pp.n{T)- The small- A: expansion 
method [2^, [1^ starts with the transformation (for 
2-D transverse plane) Pp{k) = J drpp{r) exp (ikf) , 
and examines instead the Fourier series Ppfi{k) + 
Sn>i 2pp,n(fc) cos [n(0fe — . One can then make 
a small-fc expansion of the coefficients Pp^n{k) like e.g. 
/3p,n(fc) = I]m>o'°p,™,"(^'^)™/'"'- The coefficient /Op,o(fc) 
is irrelevant to the azimuthal anisotro py and will not be 
addressed here. It is easy to see that [23[ for the n — th 
harmonic (n ^ 2), the leading term in /9p^„(fc)'s k ex- 
pansion is pp njn{— 1/2" < r"cos[n(0— 5'„)] >,<••> 
means averaging over normalized pp{f)). The case for 
n = 1 deserves special treatment: it starts with Pp,3,i 
{— 3/8 < r'^ cos (0 — ^'i) >) as one can set /Cp,i,i = by 
choosing the "right origin" of coordinates. Practically 
one then has to truncate the series, e.g. keeping only 
leading terms (assuming the relevant k is small). Due 
to the truncation, upon transforming back to the r-space 
one needs to regulate the large r part. Different from 
previous approaches, we will simply use a Gaussian fac- 
tor with a length parameter E that roughly reflects the 
scale of fluctuations. In the end, we have the following 
parameterizations for various orders of harmonics: 



PpAr) 

Pp,n(^) 



ei < r'^ > 

T S3 

n 



r \ 3 / r 



1 



! E" 



1 



7rE2 



7rE2 
, n $J 2 



,(4) 
(5) 



Here ei = - < cos {(f) - 'i' i) > / < >, e„5,2 = 
- < cos n{(j) - ■^n) > / < r" >, < r" >i/"= 
3.88, 4.19, 4.44, 4.63, 4.80, 4.95 f m for n = 1,2,3,4,5,6. 
The root-mean-radius-square V < > — 4.19 fm is de- 
fined as a and will be used as a "ruler" for E. Therefore 
by specifying the values of e„ (as well as \l/„) for an event, 
one then fixes the Pp^n 

(r) above and obtains a particular 
density profile in Eq.Q with harmonic fluctuations. 



B. The Jet Energy Loss 

In this study, we use the geometric model approach 
for computing the jet energy loss. Such models reflect 
the generic geometric features (e.g. the path-length de- 
pendence) that are most crucial for describing geometric 
data 0-11] • We assume that the flnal energy Ef oi a, jet 
with initial energy Ei after traveling an in-medium path 
P can be parameterized as Ef = Ei x fp with the sup- 
pression factor fp given by: 



fp = exp 



s{i)]s{i)rdi 



In the above the s{l) is the entropy density of local mat- 
ter at a given point on the jet path, while the k(s) is the 
local jet quenching strength which as a property of mat- 
ter should in principle depend on the local density s{l). 
There can be different choices of the parameter m for 
path-length dependence (e.g. LPM-motivated quadratic 
or AdS/CFT- motivated cubic) and of the jet-medium in- 
teraction k{s). In the present study, we use the near-Tc 
enhancement model as in Q] [1| ; assuming m — 1 and in- 
troducing a strong jet quenching component in the vicin- 
ity of Tf, (with density Sc and span of s-u,) via 



k(s) ^ K[l + ^exp{~{s - Scf/sl)] 



(7) 



with ^ — 6, Sc — 7//m3, and = 2/ frri^. (see [3| 
for details.) As aforementioned, the parameter k will be 
flxed by Raa ~ 0.18 in the - 5% collisions at RHIC 
= 200GeV. 

It shall be emphasized that the jet path P is deter- 
mined by the initial jet spot coordinates on the trans- 
verse plane as well as the the azimuthal angle 4> for the 
transverse orientation of its propagation. After averaging 
over all jet paths (including all possible start points prop- 
erly weighed by binary collision density and all equally 
possible orientations) one may then obtain the Raa- 



Raa =< (/p)"-' 



>p 



(8) 



(6) 



where the exponent n comes from measured reference p- 
p spectrum (see e.g. for a detailed account). The 
value of n depends on collision energy and throughout 
this paper we focus on the RHIC y/s — 200GeV collisions 
with n « 8.1. 

The azimuthal angle (^-dependence can be studied by 
averaging over jet paths with a particular azimuthal ori- 
entation, i.e. RaaW =< {fp)"~'^ ^P{<j>) f^o^^ which 
various harmonic components are derivable by Fourier 
decomposition. In this study for each given geometry, 
we compute Raa{4>) as a function of azimuthal angle <fi 
by integrating over all initial jet spots on the transverse 
plane weighed by the binary collision density at each spot 
for a given value of (j). 



III. JET RESPONSE TO HARMONIC 
FLUCTUATIONS 

Here we report our results for the jet response to har- 
monic fluctuations: see Fig. [TJ The results are obtained 
with E = (7 = 4.19 fm, which indicates the scale of the 
fluctuation is roughly the same as that of the isotropic 
background. In our approach, the fluctuations e„ are 
assumed to be small so that the jet response to density- 
fluctuations is approximately linear. We will flrst demon- 
strate that indeed for n — th harmonics alone with am- 
plitude e„, the jet response to a good approximation de- 
pends on e„ linearly (hence we can deflne = Vn/tn)- 
The solid curves in each panel clearly show that for 
n = 1,2,3, the relations between u„ and e„ are linear. 
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FIG. 1: (color online) The calculated vs. e„ for various harmonic fluctuations n = 1,2, 3, 4, 5, 6 with parameter E=(t. Three 
different responses are shown: the case with only pp fluctuation ( "Part" , green dashed curves) , the case with only pc fluctuation 
("CoU", blue dotted curves), and the physical ceise with both ("Tot", red solid curves) (see text for details). 



while for n = 4, 5, 6 only very mild nonlinearity starts to 
develop near e„ = and the response stays quite linear 
for the physically relevant region 0.05 < e„ < 0.15. For 
each harmonic fluctuation, we have checked that the ini- 
tial axis angle '^n of fluctuation coincides with the final 
axis angle 5';^ of the jet distribution anisotropy. We've 
also varied the parameter E and found that the linearity 
gets even better with larger S and strong nonlinearity 
develops for small E. 

One important issue is to understand the different in- 
fluences on the jet response due to the fluctuation in the 
participant density (which mostly concerns the bulk mat- 
ter that quenches the jets) and that in the binary collision 
density (which concerns the profile of initial jet spots). 
The difference could be tricky and counter-intuitive (see 
e.g. some discussions in [33|). Here we clearly demon- 
strate the difference by showing three different responses 
in Fig. [H the case with only pp fluctuation ( "Part" , green 
dashed curves) , the case with only p^ fiuctuation ( "Coll" , 
blue dotted curves), and the physical case with both 
("Tot", red solid curves). As one can see, the contribu- 
tions to the jet response due to fluctuation in participant 
density and that in collision density differ drastically and 
even have opposite signs (and therefore cancelation) for 
n — 2,4,6. It is clear that for event-by-event studies of 
jet response, both fluctuations have to be fully taken into 
account. 

Finally in the Tab. U we show the jet response coef- 
ficients Xn extracted from the physically most relevant 
region 0.05 < < 0.15 where we also have linear de- 



pendence for all harmonics. The results as a response 
spectrum vs. n are also plotted in Fig. (3] The dom- 
inant responses come from the first harmonics. We also 
see certain sensitivity to the parameter E: the smaller E 
(smaller scale for fiuctuations) leads to stronger responses 
and steeper decrease of Xn "^itti n. 



Total Response 


x1 


x5 


X§ 


x1 


xt 


xt 


E = a 


0.52 


0.15 


0.080 


-7.7e-4 


-3.5e-3 


-5.2e-4 


S = 1.25(7 


0.19 


0.13 


0.036 


6.6e-4 


-5.1e-4 


-9.5e-5 



TABLE I: Total responses to various orders of harmonic fluc- 
tuations in the initial condition using different values for the 
parameter E (see text for details). 



We end this section by discussing certain check we've 
done. One issue is about possible "mixing" or "interfer- 
ing" effect among different harmonic fluctuations when 
they all co-exist with varied respective main-axis for each. 
This should be clarified as in reality each collision event 
would automatically come with all these harmonic fluc- 
tuations. We've done the following test: we include on 
top of the isotropic background all the n — 1,2, 3, 4, 5, 6 
harmonic fluctuations simultaneously with randomly as- 
signed e„ (in the linear regime though) and randomly 
chosen axis \E'„ (different for each), and then calculate 
the jet quenching result i?^^ (</>). We've found that 1) 
the main-axis for each harmonics ^P^^ as determined by 
maximizing the corresponding cos((/) — ^';{) component 
in flnal state RAAi4') agrees with the initial from 
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FIG. 2: (color online) The jet response Xn ~ to various 

orders of harmonic fluctuations in the initial condition, with 
the blue solid curve for parameter — a and the red dashed 
curve for parameter E — 1.25a. (see text for details). 



fluctuation with less than 1% difference; 2) the response 
Vn/^n also agrees with the values we extracted above by 
treating each harmonics separately. Another issue that 
we've checked is how such response would change with 
different jet energy loss models. We've also done the 
calculation for two other models: the model with path- 
length-square dependence and constant jet-medium in- 
teraction, i.e. k{s) — K and m = 1 in Eq.([6l); and the 
model with path-length-cubic dependence and constant 
jet-medium interaction, i.e. k{s) = k and m = 2 in 
Eq.(l6]). The response coefficients Xn from different mod- 
els are not drastically different, show similar decreasing 
trend with increasing n, and may become distinguishable 
when realistic and accurate data comparison can be done. 



IV. THE "HARD-RIDGE" IN DI-HADRON 
CORRELATIONS 

In this section, we study the correlation between the 
hard and soft sector due to their respective correlations 
to the common initial condition with fluctuations. 

Suppose for an arbitrary central collision event, we 
have the initial azimuthal distribution of matter density 
fluctuations characterized by a series of harmonics with 
varied (and uncorrelated) magnitude e„ and axis orienta- 
tion ^fn. According to our results on the jet response, we 
would expect the final state high-pt hadron distribution 
to be of the form 



dN 



hard 



dydcj) 



E 

n=l,2,3,. 



2«^cos[n(0-*n)] (9) 



The bulk matter, on the other hand, will generate the 
harmonic flows during hydrodynamic expansion and lead 
to the flnal state low pt hadron distribution of the form 



dydcj) 



1 



E 

n=l,2,3, 



2<cos[n((/. (10) 



Note that for a given event both the soft and the hard 
responses to each order of the harmonic fluctuations are 
commonly aligned to the same corresponding angle 
from the initial condition — such a correlation between 
the hard and soft responses will survive the average over 
many events. 

We can therefore define a hard-soft correlation function 



C[A0] = 
/t^^2^<5(02 



01 - A0) < 



dN' 



oft 



dydtpi dyd<p2 



> 



2ir 



i(p2 dN^^ft 
It: dyd(l)2 



> 



-1(11) 



where the < > means averaging over events. Since the 
initial harmonic fiuctuations have their orientation angles 
^'n vary randomly from event to event and uncorrelated 
among each other (at least so in central collisions), we 
obtain from Eg . (flUl) the following 



E 



> cos(nA0) 



(12) 



In the linear response approximation we can use t)^'" = 
X^'^^n to further get 

C[A^]^ J2 2x''x' < > cos(nA0) (13) 

n=l,2,3,... 

We therefore see that the common correlations of the 
hard and soft responses with various harmonic fluctua- 
tions indeed lead to the hard-soft correlation that sur- 
vives the event average. 

It it tempting to quantitatively examine the fea- 
tures of such hard-soft azimuthal correlation. To 
do that, we compute C[A0] by including the n — 
1,2,3,4,5 harmonics in Eg. p^ . For the initial fluc- 
tuations we use results from Monte Carlo simulations 
of initial conditions for most central collisions, which 

are ^6^^12,3,4,5 ~ 0.05, 0.08, 0.08, 0.07, 0.07 respectively 
(see e.g. (l^jl^jfl^ji^). For the soft response coeffi- 
cients we take the results from hydrodynamic modeling, 
which are Xn=i,2^Afi ~ 0.2, 0.26, 0.21, 0.14, 0.086 respec- 
tively (see e.g. |22|[23] [25] 26]). The hard response coef- 
ficients are taken from our results: since such coefficients 
depend on the parameter S we will show results for both 
values of E as studied in the previous section. The so- 
obtained correlation function C[A(/)] is plotted in Fig. 3: 
the blue solid curve is obtained using Xn results from 
Y, ~ a case while the red dashed curve is from E = 1.25a 
case. A few remarks are in order: 

1) there is a robust and narrow near-side peak around 
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-10 12 3 4 

FIG. 3: (color online) Hard-soft correlation C[A0] from re- 
sponses to common geometric fluctuations in the initial condi- 
tion, with the blue solid curve using jet response results from 
E = (T case while the red dashed curve from E — 1.25(j case 
(see text for details). 

A(j) = 0. Actually for both cases we've fitted very well 
the peaks by a Gaussian function with width about 0.8 
radian which is very close to the experimental data. This 
main feature, we believe, provides a natural explanation 
of the "hard-ridge" structure as arising from superpo- 
sition of multiple harmonics (dominantly from the first 
three harmonics). 

2) the away-side structure around = tt is sensitive to 
the parameter S controlling the overall size of fluctua- 
tions. Depending on that parameter, the relative ratios 
among Xn change and as a result the away-side shape 
could change between a simple broad peak and a double- 
hump structure. 

While certainly interesting to do, a quantitative compari- 
son between the obtained correlation and the experimen- 
tal data is not ready yet, due to a number of approx- 
imations used in the modeling and the complication in 
experimental methods (e.g. observable definition, trigger 
and associate pt selection, the ZYAM procedure, etc). 
Some of these issues will be discussed at the end. 

One might also think about possible measurements of 
hard-hard correlations, i.e. di-hadron correlations with 
both hadrons' pt to be high enough that they are most 
likely from jets. Such measurements might become feasi- 
ble at LHC energies when enough events with more than 
one pair of jets could be collected. Along similar consid- 
eration as before we may expect a component in the hard- 
hard azimuthal correlation due to geometric fluctuations 
of the form - E„=i,2,3,... 2 ix'"? < > cos(nA0). 

A remnant of a near-side ridge might be visible, while the 
usual di-jet back-to-back correlation will be dominant. 



V. SUMMARY 

In summary, we've studied the jet response (partic- 
ularly azimuthal anisotropy) as a hard probe of the 
harmonic fluctuations in the initial condition of central 
heavy ion collisions. By implementing the fluctuations 
via cumulant expansion for various harmonics quantified 
by e„ and using the geometric model for jet energy loss, 
we've computed the response Xn — Vn/^n- Combining 
these results with known results for hydrodynamic re- 
sponse of the bulk matter expansion in the literature, 
we've shown that the hard-soft azimuthal correlation 
arising from their respective responses to the common 
geometric fluctuations reveals a robust and narrow near- 
side peak which may be the explanation of the "hard- 
ridge" observed in experimental data. 

We end by discussing a number of issues that are not 
fully addressed in the present paper and will be further 
investigated. The approach of quantifying fluctuations 
by cumulant expansion with one harmonic fluctuation at 
a time has the advantage of clearly demonstrating the 
response of jet energy loss to each order of harmonics, 
but for realistic modeling and comparison with data one 
needs to use Monte-Carlo generated initial matter pro- 
files born with various harmonics and more irregularity 
— this is currently being pursued and will be in future 
publication. Using the Monte-Carlo initial conditions 
will also allow extending our study to non-central col- 
lisions and eliminating the ambiguity of the away-side 
structure. At present the various jet energy loss models 
suffer from uncertainty for the energy loss in the pre- 
equilibrium matter Q Q (38| , and it would be interesting 
to find ways of using geometric features of jet quenching 
to constrain such uncertainty. While explanation of the 
near-side peak structure in the hard-soft di-hadron az- 
imuthal correlation by the harmonic fluctuations appears 
to be robust, the away-side structure is rather sensitive to 
the details of the fluctuations and its interpretation niay 
be tricky due to e.g. transverse expansion dynamics[39| 
and background effects like transverse momentum conser- 
vation 40;], which all require further scrutiny. All that 
said, we emphasize again that the hard probe of geo- 
metric fluctuations provides a new and useful tool for 
exploring both the initial condition and the jet energy 
loss mechanism in heavy ion collisions. 
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